rcAddPt <- function(pt,xrange=c(-5,5),yrange=c(-5,5),col=2){
segments(pt[1], yrange[1], pt[1], yrange[2], col='gray',lty=2)
segments(xrange[1], pt[2], xrange[2],pt[2], col= 'gray',lty=2)
segments(xrange[1], yrange[1]+diff(pt), xrange[2],yrange[2]+diff(pt), col= 'gray',lty=2)

segments(pt[1], yrange[1], pt[1], pt[2], col=col,lty=1,lwd=2)
segments(xrange[1], pt[2], pt[1],pt[2], col=col,lty=1,lwd=2)
segments(pt[1], pt[2], xrange[2],yrange[2]+diff(pt), col=col,lty=1,lwd=2)
points(pt[1],pt[2],pch=19,cex=1.5,col=col)
}


xrange <- c(-5,5)
yrange <- c(-5,5)

pt1 <- c(1,3)
pt2 <- c(-1,2)


plot(xrange,yrange,col='white',xlab=expression(u[1]),ylab=expression(u[2]))
rcAddPt(pt1,xrange,yrange,'red')
rcAddPt(pt2,xrange,yrange,col='darkblue')

# -----------------------------------------------------
ln <- 100
x <- seq(xrange[1],xrange[2],len=ln)
y <- seq(yrange[1],yrange[2],len=ln)